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Abstract 

The processes of lepton-nucleon scattering, including ones with both polarized 
beams, at high energy provide relevant information about interaction and particles 
structure, allowing to analyze nucleon spin structure. As energy and experimen- 
tal accuracy rise, necessity to improve Born cross sections and polarized asymme- 
tries with higher order radiative corrections becomes substantial. In this report we 
stress on lowest order bremsstrahlung corrections treatment using helicity ampli- 
tudes method as applied to actual nowadays charged current lepton-nucleon deep 
inelastic processes, that allows to simplify matrix element calculation procedure. 
Real photon emission contribution is calculated by means of Lorentz-invariant for- 
malism. Kinematical peculiarities on bremsstrahlung correction are discussed. 

1 Introduction 

The processes of deep inelastic lepton-nucleon scattering (DIS) are of interest at 
present and planned experiments, today with particular emphasis on both polar- 
ized beams interaction investigation, as it provides essential data on the internal 
structure of the nucleon spin. Special interest to charged current interaction is 
connected with the absence of large electromagnetic effects contribution to these 
processes. Some extensive reviews on nowadays and forthcoming experimental fa- 
cilities on such processes can be found for instance in refs. [TJ [2]. Asymmetries 
withdrew from phenomenological interaction parameters on certain experiments al- 
low to extract detailed information on nucleon's spin, concealed in polarized struc- 
ture functions g± and 52(5,6) or individual quark contributions to nucleon's spin. As 
expected, obtained information can be used to expand and to refine nucleon nature 
knowledge, to compare experimental data with other related experiments on nu- 
cleon structure (e.g. neutral current or pure electromagnetic DIS ones, which have 
been studied at a stretch of many years, see for instance refs. [3-5]) as well as with 
Standard Model predictions or perhaps to search deviations from it. 

Processes in question have been investigated before mainly at Born approach, for 
instance in our previous papers (see in refs. (6JE]) we realized Born level phenomeno- 
logical analysis in comparison with quark-parton model approach, Born asymmetry 
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analysis with stress on polarized structure functions extraction scheme. In this 
report we restrict oneself to detailed treatment of the bremsstrahlung correction 
calculation, as correct treatment with observed experimental data at high energies 
requires allowance for various radiative effects. Here to perform calculations we use 
the formalism of helicity amplitudes method offered firstly in ref. [HJ [9] relevant 
for single and multiply bremsstrahlung processes. Using of such analytical method 
allows to practically avoid intermediate operations with traces of Dirac matrices 
products and undesirable calculations of cross-elements of S'-matrix. This method 
of matrix element calculation mainly consists in special representation of 4-vectors 
of the photon polarization through expressions with bispinors, in using of special 
transformation rules likewise Chisholm identities and in treatment with u^(p)u±(k) 
constructions as simple scalar function of p and k to cancellate unnecessary terms. 



2 Radiative corrections 

To calculate radiative corrections we employ quark-parton model, which allows to 
obtain reasonable quantitative predictions for nucleon and leptonic bremsstrahlung 
contributions. Feynman diagrams of processes in question 

1(1) + N — v{p) +X, 
v{v)+N -+l(T)+X, 

are presented in FIG.l (one of the diagram in particular case vanishes, as neutrino 
contain no charge). 
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Figure 1: Photon bremsstrahlung diagrams for charged current lepton-nucleon DIS pro- 
cesses. 

Let's consider firstly the case of l{ = I, If = vi, qi, qj. To make use of helicity 
amplitudes method one should represent photon polarization vectors as following: 



e± = f± = N q 



q"q'k{l T 75) " kq"q'{\ ± 75 ) 



W-(q"q'Wk)(q>k) 



e ± = 2N a 



q"q'k(lTK)-kq"q , (l±K)T 
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where q • ,q ' ,PpPf correspond to incoming I or outgoing v\ leptons and quarks qtqj 
momenta, k is emitted photon momentum. This representation conserves common 
requirements e^e^ = 0, e^e^* = e^e^* = 1, e^k = 0. Given form suitable for 
calculation of the leptonic bremsstrahlung term. For hadronic terms one should 
involve form, dependent on quark's momenta as following 



Op 



p'}plk(lTj5)-kp'fPi(l±J 5 ) 

{p" s pWfk){my\ 



+ 



where 



o±i<P 



N p N q Sp 
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p" f p' i kq"q'k{l =F 75 ) 



Free parameter (3 can be omitted, as longitudinal component. 

One can get the following expressions for matrix element using technique thor- 
oughly described in ref. [9] 



M_ 



AN p e +i ^e f 



u(q")q'<P^(Pf)Pi x u(q')D w , 



2ky l0 p» Q 

A[N q + N p e~^e, 



i2ky ioP " fQ 

xu(p'})q"u(q')D W , \A\ 2 



-u(q")p'f u (p'i) x 



e 2 G 2 F M§ v 
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2 HQo 

Here D w - W-b oson propagator, signs ± refer to the particle helicity and photon 
polarization in the following order: /, f/, 7, qi, qf, and ej - initial and final quark's 
charges. 

The advantage of using helicity amplitudes method is that it allows to obtain 
squared matrix element directly without interference terms common in straight- 
forward calculations. It is readily to show, that squared matrix element of the 
lepton-quark process Iqi —* viqf with real photon emission have the following form 



|M_ 



+ |M_ 



Ae 2 G 2 F 



q qokoP i0 Pfo 



{Q 2 + M^Y 



ei(q'k)(q"k) 



( 4 


p'fX 


{P'lk 


p»k)_ 



q'k q"k 
n„n\2 



+ 



{p" f q" 



ip'ik)ip" f k){q'k){q"k) 

e f (q'q")(pW) 2 \ 
(pJfc)(P/fc) J 
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Similar formulas for other cases of electroweak lepton-quark scattering can be 
evolved using transformation rules 



Iqi -> viqf 
Iqi -> &iqf 
Iqi -> P/<?/ 



l'<->q", p'i^p'f, e»<->e/. 



We use here the following common notations for kinematical variables: 
Ql = q'-q", Qf^-2q'q", 

Q h = Qi-k = p" f - qi « -2 P y;, 

5, = -2^, 

u = -2p[k, z 2 = -2q"k, 
v = -2p" f k = u - Q 2 + Ql, 
Zl = -2q'k = z 2 - Qj + Q\, 

where x^q = —Ql^/2p'Qh, Uh[i] = ~^p'Qh[i}/ S - standard hadron and lepton scal- 
ing variables, p' and p" - nucleon and jet 4-momenta. To obtain cross sections or 
polarized asymmetry including radiative corrections one should switch from lepton- 
quark interaction to lepton-nucleon one integrating over quark momenta being car- 
ried in the nucleon, and over radiated photon momentum. If we suppose quark to 
possess momentum p[ = x^p' with the probability of f(xih), the first integration 
over p\ could be performed by means of the following substitutions: 

f{Xlh) J 2^Q~ h = Si XlhS XhS = ^' 

X % - x lh X -^ Xh X = ^(l- y h )S, 

Vh 

u^x ih u^x h u, v -> x h u - Qi + Ql, 

keeping Ql = (Qi — k) 2 , Qf, z\ and z 2 unaltered. Here S = —{jpl + q') 2 . 

To integrate over photon momentum one can use covariant method of integration 
described, for instance, in ref. [101111]. permitting to integrate directly over Lorentz- 
invariant kinematical variables. Covariant calculation has advantage of missing of 
the sophisticated Monte-Carlo techniques but presence of the analytical integration 
as well as it can be carried out for various kinematical experimental configurations. 

Firstly lets imply the following suitable phase space transformation, allowing to 
derive from its general form the expression containing introduced before invariant 
variables: 



dM 2 dQl^^6[{ Ql -k) 2 + M 2 



h 



, ,„ 9 , ,„ 2 dz 



:<5 [Ql ~ [P"-P'f] = —d yi dQ 2 dy h dQ 



2 ah ~* h 2jR z 



Here R z is the Gram determinant [12] of 4- vectors q', p', q", p" 

R z = -A i (q',p',q",p") = 
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which can be expressed as quadratic polynomial of z\ or z 2 variables defined before 

R z = -Az 2 + 2Bz - C, 
where the coefficients in the ultrarelativistic limit are 

Ai,2 = y 2 s 2 + AM 2 Ql 
B 1 = -2M 2 Qf{Q 2 - Ql) + ( yi Q 2 - y h Q 2 )S 2 + 
+(1 - yi)S 2 QUvi ~ Vh) - m 2 {2M 2 Q 2 h + 2M 2 Q 2 - S 2 y hVl ), 
B 2 = 2M 2 Q 2 (Q 2 - Q 2 ) + (1 - yi )( yi Ql - y h Q 2 l )S 2 + 
+S 2 QHvi ~ Vh) - m 2 (2M 2 Q 2 + 2M 2 Q 2 - S 2 y hVl ), 
C 1 = S 2 [Q 2 h + (-l + y l -y h )Q 2 } 2 + 

+Am 2 Q 2 ( yi -y h f{l-yi), 
C 2 = S 2 [(l-y l )Q 2 h -{l-y h )Q 2 l } 2 + 
+Am 2 Q 2 {y l -y h f{l-y l )-\ 

In presented above expression we simplified common phase space by means of aux- 
iliary invariant variables z\ or z 2 . Next one can employ the following integration 
scheme: 

max ^fimax ^max „ 

dz dyi dQf 



da ~ J dy h J dQ\ J 



Vhl 



or 



? 2 

* h min 



!/l max ^ I max 2 m ax _ 

Vl min O? - ^min 

^ / mm 



y h Sy/R~ z 



-A 



A = \M\ 2 fi(x h ,Q 



and so on, dependently on desired final variables. Here matrix element \M\ ex- 
pressed in terms of Q 2 ! h uyi[h] and Zip] have the following form: 



\M(Si,Xi,Q h ,Q h z [li2] )f 



e 2 S 2 Q 2 2 ; 

+ ( Z2 - Ql - Xl ) 



[Q 2 h uv + (e 2 Qfziz 2 - eiQ 2 u(zi + z 2 ) + eiQ 2 (u - v)(SiZ 2 - X^))] 



2uvz\z 2 



for k = 1, If = u x , qi, qf and 



{MiSuX^Q^Q^z^f 



Ql x f , e}Qi X ? , e 2 Q 2 {z 1 + Q 2 -S % f 




+ 



e f Xf [uQfjzi + z 2 ) - (u - v){SjZ 2 - XjZi)] 



2uVZ\Z 2 



for k =1, If = u h q h q f . 

Calculation of the integral over z can be carried out using simple table integrals 
of the form 



To evaluate remaining integrals one should firstly use some parameterizations on 
quark distribution functions fi(xh,Ql), for instance QCD-based fits from ref. [13], 
and then choose final variables in which result cross-section or asymmetry will be 
expressed. 

In order to calculate infrared contribution, one can apply the limits lim (u — v) = 
and lim (£2 — Z\) = to presented above relations. When calculating this part 



requires using of some regularization method as being infrared divergent part, for 
instance by introducing virtual photon mass or applying dimensional regularization 
method. Here we didn't stress on soft photon emission contribution calculation, 
as virtual loop contribution and soft emission one, infrared divergent separately, 
compensate partly each other, remaining uncompensated part have minor influence 
on the asymmetries for processes in question as reduced factor. 

Kinematical peculiarities and variables limits for these integrals are thoroughly 
described e.g. in ref. [12], here we give only kinematical relations, necessary for 
determination of the integration bounds. Imposed constrains on the physical region 
of invariants have the following form in terms of kinematical A- functions (see [12] ) 



\{\sAu\) < 0, \{\ q ,\ h ,\ k ) < 0, 

X(X s ,X T ,X h ) <0, 
A fc = (yi - y h ?S 2 , \ S = S 2 - 4m 2 M 2 , 
X^il-ytfS 2 , X q = yfS 2 + 4M 2 Q 2 , 
X h = y 2 h S 2 + AM 2 Q 2 , X T = (1 - y h ) 2 S 2 - 4M 2 z 2 , 



The boundary equations for this conditions can be expressed in the form of three 
equations 



consequently. 

One can obtain certain integration bounds for chosen final variables by combining 

,1 . , -, n min,max j min,max ■ r j-,- 7-. ^ n 

these constraints with z 1 and z 2 emerging from the condition H Zl 2 > 0. 

In ref. [6l [7] we numerically calculated hard real photon emission contribution 
with low photon energy cut parameter e cu t, where unpolarized quark distribution 




n 



= -2... 2. 



X(x, y, z) = x 2 + y 2 + z 2 — 2xy — 2yz — 2xz. 



M 2 (Q 2 - m 2 ) 2 + Q 2 Vl S 2 - m 2 S 2 yi (l - m ) - Q 2 X S = 0, 

y 2 S 2 Q 2 h + y 2 S 2 Q 2 - M 2 (Q 2 - Q 2 ) 2 - Vl y h {Q 2 + Q 2 )S 2 = 
(2M 2 z\ - 2m 2 M 2 + y h S 2 + 2M 2 Q 2 h ) 2 - X s X h = 0, 
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Figure 2: Born (solid line) and next to Born order (dotted line) longitudinal polarized 
asymmetry A\\(xi,y{) for IN —>■ viX at E = 100 GeV and e cnt = 1 MeV. 



functions from ref. [13] with polarized one from ref. [H] were taken. In FIG. 2 we 
cite comparison of the corrected polarized asymmetry 

A _ U 2 a^ fgW\ ffoM | rfVn 
" \ dxdy dxdy ) \ dxdy dxdy J ' 

d?a/dxdy - differential cross section, f (J,) corresponds to lepton helicity value 
—1(1), ff (JJ-) for nucleon spin, parallel (antiparallel) to lepton momentum, with 
the Born one in lepton scaling variables x\ and y\. As for relative contribution, it 
increases with x\ growth, reaching values for xi = 0.9 approximately larger by 20 % 
than for X[ = 0.01; it roughly inversely exponentially dependent on yi, not exceeding 
~ 30% for yi > 0.3 and x { < 0.9, ps 20% for y t > 0.3 and x t < 0.25 and w 10% for 
y; > 0.3 and x; < 0.1. Asymmetry absolute values depend on y\ noticeable unlike 
Born ones. 

As for multiply bremsstrahlung, such contributions can be obtained on the ba- 
sis of single one by applying renormalization group equations, and will be treated 
hereafter. 

In conclusion, real photon bremsstrahlung contribution calculated here affects 
significantly cross-sections and polarized asymmetries, so it's necessary to take into 
account such contribution, in that way it can be used for certain future experimental 
needs with the aim of precision DIS data extraction to improve accuracy of quark 
distribution functions as well as to detail nucleon's spin structure. 
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